Propagation Failure in Excitable Media 
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We study a mechanism of pulse propagation failure in excitable media where stable traveling 
pulse solutions appear via a subcritical pitchfork bifurcation. The bifurcation plays a key role in 
that mechanism. Small perturbations, externally applied or from internal instabilities, may cause 
pulse propagation failure (wave breakup) provided the system is close enough to the bifurcation 
point. We derive relations showing how the pitchfork bifurcation is unfolded by weak curvature 
or advective field perturbations and use them to demonstrate wave breakup. We suggest that the 
recent observations of wave breakup in the Belousov-Zhabotinsky reaction induced either by an 
electric field Q or a transverse instability are manifestations of this mechanism. 



I. INTRODUCTION 

Failure of wave propagation in excitable media very 
often leads to the onset of spatio-temporal disorder. In 
the context of electrophysiology it may lead to ventricu- 
lar fibrillation. Numerous studies have appeared in the 
past few years demonstrating conditions and mechanisms 
for failure of propagation in excitable and bistable me- 
dia |^-^| • Failure may occur by external perturbations or 
spontaneously by intrinsic instabilities. Experimental ex- 
amples include wave breakups in the excitable Belousov- 
Zhabotinsky (BZ) reaction induced by an electric field [Q 
and by a transverse instability [ 2j . 

Recently we have attributed domain breakup phenom- 
ena in bistable media to the proximity to a pitchfork front 
bifurcation illustrated schematically in Fig. 0a ||To|jLl|| . 
As shown in Fig. |l|b, near the bifurcation small per- 
turbations, like an advective field or curvature, unfold 
the pitchfork bifurcation to an S'-shaped relation. If the 
perturbation is large enough to drive the system past 
the endpoint of a given front solution branch the front 
reverses direction. A local reversal event along an ex- 
tended front line in a two-dimensional system involves 
the nucleation of a pair of spiral waves and is usually 
followed by domain breakup. The front bifurcation illus- 
trated in Fig. ^]a has been referred to in the literature as 
a Nonequilibrium-Ising-Bloch (NIB) bifurcation |l^ , p^[| . 

In this paper we extend these ideas to wave breakups in 
excitable media. The NIB bifurcation is replaced in this 
case by a subcritical pitchfork pulse bifurcation as shown 
in Fig. ||a. The typical unfolding of that bifurcation is 
shown in Fig. ^b. Similar to the case of front solutions 
in bistable systems, perturbations that drive the system 
beyond the edge point of a given pulse branch may ei- 
ther reverse the direction of pulse propagation, or lead 



to pulse collapse and convergence to the stable uniform 
quiescent state (not shown in Fig. ||b). The convergence 
to a uniform attractor is more likely to occur for pulse 
structures than for fronts, and has always been observed 
in our simulations. Thus, in our study, the critical value 
of the bifurcation parameter, e/, at which the upper and 
lower branches in Fig. ||a terminate, designates failure 
of propagation. Reversals in the direction of propaga- 
tion (rather than collapse) have been observed in exper- 
iments on the Belousov-Zhabotinsky reaction subjected 
to an electric field fl4|| . 

In two space dimensions a local collapse of a spatially 
extended pulse amounts to a wave breakup. By numeri- 
cally integrating a FitzHugh-Nagumo (FHN) type model, 
we demonstrate two scenarios of wave breakups: breakup 
induced by an advective field, modeling an electric field in 
the BZ reaction, and breakup induced by a transverse (or 
lateral) instability. The two scenarios, observed both in 
experiments [0j|] and numerical simulations , reflect 
the same mechanism: the ability of weak perturbations 
to drive transitions from one of the pulse branches to the 
uniform attractor when the system is close to failure of 
propagation. 

In Section II we describe the derivation of a pulse bi- 
furcation diagram for an FHN model. The derivation ap- 
plies to both excitable and bistable media. The informa- 
tion contained in this diagram is used to draw the propa- 
gation failure line in an appropriate parameter space. In 
Section III we consider the unfolding of the pulse bifurca- 
tion by an advective field and demonstrate wave breakup 
near the propagation failure. The unfolding by curvature 
is studied in Section IV and wave breakup induced by a 
transverse instability is demonstrated. We conclude in 
Section V with a discussion. 
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FIG. 1. Bifurcation diagrams for fronts in the bistable 
FHN reaction-diffusion system, a) At the Nonequilibrium 
Ising-Bloch bifurcation a stationary front becomes unstable 
to a pair of counterpropagating fronts as a control parameter 
e is varied. The solid lines represent a branch of front so- 
lutions with speed c. b) Unfolding the bifurcation near the 
critical point, e c , gives an S-shaped relation in the unfolding 
parameter P. 
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FIG. 2. Bifurcation diagrams for pulse solutions in the ex- 
citable FHN system, a) A typical relation for the pulse speed, 
c, vs the system parameter e gives a subcritical pitchfork bi- 
furcation, b) Unfolding the bifurcation near the critical point, 
e c , gives a multiple S-shaped curve for c in the unfolding pa- 
rameter P. 
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II. A BIFURCATION DIAGRAM FOR PULSE 
SOLUTIONS 

We derive the pulse bifurcation diagram using an 
activator-inhibitor model of the FHN type. We assume 
that the activator varies on a time scale much shorter 
than that of the inhibitor and use singular perturbation 
theory ||l5|-|l7|] . Specifically, we study the pair of equa- 
tions 



u t = e 1 (u-u 3 — v)+5 1 V 2 u, 
v t = u — a%v — ao + V 2 v , 



(1) 



where [i = e/S <C 1 and the subscripts t denote partial 
derivatives with respect to time. We consider a periodic 
wavetrain of planar (uniform along one dimension) pulses 
traveling at constant speed c in the x direction. Each 
of the excited (or "up state") domains occupy a length 
of A+. The recovery (or "down state") domains are of 
length A_. In regions where u varies on a scale of order 
unity Eqns. ([!]) reduce to 



u xx + cv x + u +( v ) ~~ a i v ~ °o = , 
°xx + cv x + u - ( v ) ~ a i w ™ °o = , 



-A_ < x < 0, 
0< X <A+, (2) 



where \ = x — ct and u±(v) are the outer solution 
branches of the cubic equation u — it 3 — v = 0. For 
a\ sufficiently large we may linearize the branches u±(v) 
around v — 

u±(v)&±l-v/2. (3) 
We solve Eqns. (0) using the boundary conditions 



u(-A_) = v b , 
v(0) = v f , 
v(X+) = v b , 



(4) 



where Vf and v b are yet undetermined, and the linear 
approximation ([|). The solutions are 



v — A + exp(crix) + B + exp(cr 2 x) + v + 
v — A- exp(crix) + B- exp((T 2 x) + w_ 

with 



-A_ < x < 0, 
< x < A+ , 



01.2 



A+ = 



2 

{Vb 



± 



+ ai + 1/2, 



v±) - (vf - v±) exp(=Fcr 2 A T ) 



exp(=FcriA T ) - exp(=F<7 2 A T ) 
-(v b - v±) + (vf - v±) exp(TQ"iA T ) 
exp(=FcriA T ) - exp(=Fcr 2 A T ) 



and v± — (±1 — ao)/(ai + l/2). Matching the derivatives 
of the solutions at x = 0, v'(0~) = v'(0 + ), and impos- 
ing periodicity on the derivatives, A_) = i/(A + ), we 
obtain two conditions 



<JiA+ + <T 2 i?+ = <j\A- + cr 2 i?_ , (5a) 
o\A+ exp(-CTiA_) + a 2 B + cxp(-cr 2 A_) = 

a\A^ exp(<TiA + ) + <7 2 i?_ exp(cr 2 A+) . (5b) 

Two more conditions are obtained by studying the 
"front" and the "back", that is, the leading and trail- 
ing border regions between the excited and recovery do- 
mains. Stretching the spatial coordinate according to 
C = x/\/Ji gives the nonlinear eigenvalue problems 



U CC + ai u C +u — u — Vf = 0, 



(G) 



«(0 

u(0 



u-{v s ) 
u+{v s ) 



as 
as 



-co , 



(7) 



for the narrow front region, and 

Uqq + CTjU^ + u — u 3 — v b = , 

u(C) — ► u + (vb) as C^oo, 
— * (ut) as C —>■ — oo , 



for the narrow back region. Here, 77 = \/eS. Solutions of 
(|) and (0) yield 



cr\ 



cr, = -j=v b . 



(8) 
(9) 



respectively. 

Substituting the relations (§) and (||) into Eqns. (||) 
leaves two equations for the three unknowns A+, A_, and 
c. The one parameter family of solutions describe peri- 
odic wavetrains of pulses traveling with speed 



c = C(A; 7], a , ai) 



(10) 



where A = A + + A_ is the varying wavelength of the 
family. The graph of c versus A provides the dispersion 
relation curve obtained in earlier studies |ll|[l6|]. Our 
interest here is with the behavior of a single pulse and 
therefore only the limit of large A will be considered. We 
remind the reader that in deriving these equations we 
have assumed (J, = e/S <C 1 which excludes very small 5 
values. 

We solved Eqns. (||) numerically for both excitable and 
bistable systems at large values of the period A. The 
solutions were computed by numerical continuation of 
known solutions when clq = 0, and A + = A_. They 
yield the typical bifurcation diagram for the speed c in 
terms of the parameter e as shown in Fig. ||a. At some 
critical value, e/, the branch of solutions terminates and 
past that point pulses fail to propagate. The value of e/ 
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FIG. 3. The pulse failure propagation boundary in the e — S 
parameter plane for the excitable FHN model. To the right 
of the line pulses fail to propagate. Parameters: a\ — 1.25, 
a = -0.2. 

depends on the other system parameters as well. Figure || 
shows a graph of e = e/ (5) for an excitable system. 

Note the inherent subcritical nature of the bifurca- 
tion [l8|Jl9[ |. The bifurcation becomes supercritical only 
in the limit ao — > (pertaining to a bistable medium) 
where the pulse size tends to infinity, fn that limit 
Eqns. (||) can easily be solved. The solutions are c = 
and c = ±-^y^| — r\ 2 for T) < r) c , and coincide with the 
Nonequilibrium Ising-Bloch (NIB) bifurcation for front 
solutions pfj. 
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FIG. 4. Solutions of the speed, c, vs advection, J, rela- 
tion (^|). (a) Away from the bifurcation point, variations in 
J have little effect on the pulse speed (e = 0.01). (b) Near 
the bifurcation point, small variations in J may drive the sys- 
tem past the endpoint of the solution branch and cause pulse 
collapse (e = 0.05). Other parameters: ai = 1.25, ao = —0.2, 
(5=1.0. 



III. WAVE BREAKUP BY AN ADVECTIVE 
FIELD 

Application of an electric field to a chemical reaction 
involving molecular and ionic species, like the BZ reac- 
tion, results in a differential advection pi] . Differential 
advection in the FHN equations can be modeled (with- 
out loss of generality) by adding an advective term to the 
inhibitor equation 

u t = e- 1 {u-u 3 -v) + 5- 1 V 2 u, (11) 
Vt = u — a\v — ao + J • Vv + V 2 w , 

where J is a constant vector. Looking for planar solutions 
propagating at constant speeds in the J direction, and 
rescaling e and 5 by the factor -^-j we find the influence 
on the pulse speed by the advection 

c+J = C(X; — —r],ao,ai), (12) 
c + J 

where C is defined in Eqn. (|i"o|). Figures [| show the de- 
pendence of the pulse speed c on the advection constant 



J obtained by numerically solving Eqn. (|12|) for large 
A. Away from failure of propagation (e is significantly 
smaller than e/) small variations of J have little effect 
on the pulse motion (Fig. ||a). However, close to failure, 
such variations can induce wave breakup by driving the 
system past the end point J = J c (Fig. ||b) . 

Figure || shows a numerical simulation of Eqns. ( |i"l| ) 
with J = Jx and an initial condition of a curved pulse. 
Along the pulse the effective advection field is the projec- 
tion of J onto the direction of propagation at that point. 
With J = the curved pulse propagates uniformly out- 
ward in a circular ring. Choosing J slightly greater than 
J c , a wave breakup results: the part of the pulse propa- 
gating in the x direction fails to propagate. Those parts 
propagating in significantly different directions still prop- 
agate. These results explain earlier observations of wave 
breakup induced by electric fields |jj . 
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FIG. 5. Numerical solution of Eqns. (|Ll|) with a weak ad- 

vective field J = Jx. The thick and thin lines pertain to 
u = and v — contour lines, respectively. The initial circu- 
lar pulse fails to propagate along the direction of the advective 
field and the pulse breaks. The pulse continues to propagate 
in directions different from the advective field. The equation 
parameters are the same as in Fig. ^[ 



IV. WAVE BREAKUP INDUCED BY A 
TRANSVERSE INSTABILITY 

Spatially extended pulses, like stripes or disks, may be 
unstable to transverse perturbations along the pulse line. 
Ohta, Mimura, and Kobayashi studied the case of defor- 
mations of planar and disk-shaped stationary patterns in 
a piecewise linear FitzHugh-Nagumo model p5] . Kessler 
and Levine derived conditions for the transverse instabil- 
ity of traveling stripes in a piecewise linear version of the 
Oregonator p3| . The curvature induced by a transverse 
instability can lead to the formation of labyrinthine pat- 
terns p 24 25|] or cause spontaneous breakup of a pulse 
as we will now show. 

For the model of Eqns. (3), the effect of curvature k 
on pulse propagation can be obtained from Eqn. (10) 



by rewriting the equations in a frame moving with the 
pulse pq| . Assuming the radius of curvature is much 
larger than the pulse width, and a negligible dependence 
of u and v on arclcngth and time (in the moving frame) 
we obtain 



(c + S 1 n)u' + e 1 (u-u 3 -v) 
v + (c + k)v + u — a\v — ao 



(13) 



where the prime denotes differentiation with respect to 
a coordinate normal to the front line. Rescaling e and 6 
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FIG. 6. Solutions to the speed vs curvature relation (|l4|). 
(a) Away from the bifurcation point (top), the speed c varies 
approximately linearly with the curvature k (e = 0.003). (b) 
Near the bifurcation point, small curvature variations may 
drive the system past the endpoint of the solution branch 
and cause pulse collapse (e = 0.022). Other parameters: 
oi = 1.25, a = -0.2, 8 = 2.5. 



C + S K 

c + k = C(A; r], ao, a\) . 

c + K 



(14) 



Figure ^ shows numerical solutions of Eqn. (Q) for c in 
terms of k. Far away from failure of propagation (Fig. |6^,) 
we find the usual approximate linear c — n relations for 
right (c > 0) and left (c < 0) propagating pulses ]l5| , [T7t . 
Close to failure (Fig. |6|b), small realizable curvature vari- 
ations may cause collapse. 

Equation (^) contains information also about the 
transverse stability of a pulse line. A positive slope of 
a c — k relation at k = indicates an instability of a 
planar pulse. Fig. ^ shows a simulation of Eqns. (Q) 
at parameter values pertaining to the c — k relation in 
Fig. ^b. Starting with an near planar pulse, dents grow 
due to a transverse instability. The negative curvature 
that develops induces a wave breakup. 



by the factor C+ C S +K K Eqn. (10) becomes 
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FIG. 7. Breakup of a pulse by transverse instability. The 

thick and thin lines pertain to u = and v = contour lines, 

respectively. The initial almost planar pulse is unstable to 

transverse perturbations and forms a dent. The dent grows 

and the pulse breaks at the region of high curvature. The 

equation parameters are the same as in Fig. g. 



V. CONCLUSION 

We have identified a mechanism for breakup of waves 
in an excitable media. The key ingredient of this mech- 
anism is the proximity to a subcritical pitchfork pulse 
bifurcation (as shown in Fig. ^a). Near the bifurcation 
small perturbations become significant and may induce 
failure of propagation. The nature of the perturbation is 
of secondary importance. As illustrated in Figs. f|b and 
^|b, the effects of an advective field and curvature are sim- 
ilar; they both induce wave breakup by driving the sys- 
tem past the end points of propagating pulse branches. A 
perturbation inducing breakup can be externally applied, 
like an electric field in the BZ reaction, or spontaneously 
formed, like curvature growth by a transverse instabil- 
ity. An interesting question not resolved in this study is 
the observed preference of propagation failure or collapse 
rather than reversal in the direction of propagation. 
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